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Abstract. We study the Kondo lattice model with the Hcisenberg-type RKKY- 
exchange coupling among localized /-spins in the presence of a magnetic field. By 
means of an extended dynamical mean field theory combined with the non-crossing 
approximation, we investigate the one-particle spectral function and the dynamical 
spin correlation function in the Kondo insulating phase. It is shown that the 
magnetic field and the RKKY exchange interaction both cause the instability to the 
antiferromagnetic order with enhanced transverse spin fluctuations, which give rise 
to a strong renormalization of quasi-particles as the system approaches the quantum 
critical point. This leads to a tendency to retain the Kondo insulating gap up to rather 
large fields. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.30.+h. 
1. Introduction 

Electron correlations caused by the interactions between conduction electrons and 
localized magnetic moments have provided a variety of important subjects in condensed 
matter physics. In dilute magnetic alloys, a local moment is screened by conduction 
electrons via the Kondo effect below its characteristic temperature, and the system is 
described as a local Fermi-liquid [1]. In concentrated systems, such as heavy-fermion 
systems, the indirect exchange among the localized moments mediated by the RKKY 
interaction becomes important [2,3]. The interplay between the Kondo effect and the 
RKKY interaction gives rise to competing phases in heavy fermion systems. 

Recent experimental observations of quantum critical behavior in heavy fermion 
materials have attracted renewed interest in this field. In CeCue-xAux [4] and YbRh2Si2 
[5], the system is driven to the antiferromagnetic phase by alloying, applying pressure, 
etc. In the vicinity of the quantum critical point, both Kondo and RKKY interactions 
become essential and physical properties observed in this region differ from those of 
ordinary Fermi-liquids [6] . Experimental observations of such non- Fermi-liquid behavior 
[7] stimulated detailed theoretical study beyond the standard framework [8-11] of the 
quantum critical point [12]. 
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The interplay between Kondo and RKKY interactions have been investigated using 
various methods such as slave-particle mean field theory [13,14], large- iV expansion [15], 
quantum Monte Carlo (QMC) simulations [16], dynamical mean field theory (DMFT) 
[17-20], etc. Recently an extension of the dynamical mean field theory has been 
developed [21-24], which allows us to incorporate spatially extended spin correlations. 
Si et al. have addressed the problem of non-Fermi-liquid behavior near the quantum 
critical point by the extended-DMFT (EDMFT) analysis of the Kondo lattice model 
(KLM) with the Ising-type RKKY exchange among locahzed moments [25,26]. These 
investigations focused on magnetic properties around the quantum critical point in 
heavy- fermion systems. Also, P. Sun and G. Kotliar have addressed the same issues 
by exploiting the periodic Anderson model (PAM) [27,28]. 

Magnetic fields should play a particularly important role if the system is in the 
proximity of magnetic instability, because the applied fields may possibly trigger a 
magnetic phase transition. Recent studies on the half-filled KLM in two-dimensions 
based on mean field theory and QMC simulations have indeed found that the magnetic 
field induces a second-order phase transition from the paramagnetic to the transverse 
antiferromagnetic phase [29, 30] . Also, in our previous paper we have investigated the 
field- induced phase transition in the three-dimensional case [31]. We have studied the 
half-filled PAM by means of DMFT combined with QMC. However, DMFT cannot 
incorporate inter-site fiuctuations due to the RKKY-intcraction. It is desirable to 
investigate the effects of the RKKY interaction on electronic properties by means of 
EDMFT. 

In this paper, we study the KLM with the Heisenberg-type RKKY-exchange among 
localized /-spins. By means of EDMFT combined with the non-crossing approximation 
(NCA) [32-35], the one-particle spectral function and the dynamical spin susceptibility 
are calculated. We also investigate the field-induced antiferromagnetic instabihty and 
discuss the remarkable behavior in spectral properties in the vicinity of the critical point. 

The paper is organized as follows. In the next section we briefiy summarize our 
model and method. In section. 3, we present our numerical results, and discuss the 
effects of the RKKY interaction in a magnetic field. A brief summary is given in section 
4. 

2. Model and Method 

2.1. Kondo lattice model 

We study the following Kondo lattice model including the Heisenberg-type RKKY 
interaction, 

- tY^ clcj^ -fiYl ^^i'T^^" + JK^^r Sf 

+ ^ E ■ s| - 9i^bB J2 \s-. + 5.4] , (1) 
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where Cjo- (cj^) annihilates (creates) a conduction electron on the ith site with spin o", Sd 
(S/i) is the spin operator for a conduction electron (spin-i localized /-electron). Here, t 
is the nearest-neighbor hopping matrix in the conduction band (referred to as the c band 
in the following), Jk describes the on-site Kondo coupling, and J denotes the RKKY- 
exchange among the nearest-neighbor /-spins. We consider the three-dimensional simple 
cubic lattice with dispersion — — 2t(cos A;j;-|-cos ky-\-cos k^) and the RKKY interaction 
Jq = J{cosqx + cosqy + cosqz). We focus on the Kondo insulating phase at half-fiUing, 
and take half the band width as the energy unit. The magnetic field B apphed along 
the z direction is coupled to both of the c and / electrons. We assume that the g factors 
of the c and / electrons are the same, and set g/iB — 1- We use the Kondo coupling 
Jft: = 0.5 in following discussions. 



2.2. Extended- DMFT for the Kondo lattice model 

DMFT is a powerful framework to study strongly correlated electron systems. Since its 
power was proved by providing the first unified scenario for the longstanding probrem of 
the Mott transition in the Hubbard model [36] , various correlated electron models have 
been studied by using DMFT [17-20,34,36-45]. EDMFT is an extension of DMFT so 
as to allow dynamical treatment of inter-site interactions, which can be dealt with only 
at the mean field level within the standard DMFT framework. EDMFT has successfully 
applied to the Kondo lattice model [26,46,47], the periodic Anderson model [27], the 
extended Hubbard model [48], and the t-J model [35], etc. 

In EDMFT the original lattice model is mapped onto an effective impurity model 
with a bosonic bath in addition to a fermionic bath that is included in the standard 
DMFT. The bosonic effective bath describes a fluctuating magnetic fleld and allows 
dynamical treatment of inter-site spin correlations. This approximation becomes exact 
in the limit of inflnite dimensions, d oo, provided that t and J are scaled as t/\^ 
and J/\/d respectively. It may give a good approximation even in three dimensions. 
We assume the paramagnetic Kondo insulating phase without long-range order. In this 
case, the Kondo lattice model (1) can be mapped onto the following Bose-Fermi Kondo 
Hamiltonian [23] as an effective impurity model: 

-f^EDMFT = -f^loc + -f^bath + -fornix, (2) 

Hioc - + JkS' ■ Sf - BSl - BSf, (3) 

a 

-f^bath = ^ -E^feCToLajtCT + ^(^qa^ga^a, (4) 

-fornix = V ^ ^C^aka + Ofcj^Cf^j + 

k,(7 q,a 

Here, ifioc is the local part of the effective Hamiltonian at the site chosen as an 'impurity', 
which includes the Hartree contribution of the RKKY-interaction among the /-spins, 
B = B — 2Jq=[Qfi^Q]{S^)S^ . The effective baths are described as i/bath) where a^^ creates 
a fermionic bath with the dispersion E^a, and hga (a = x,y,z) creates a bosonic bath 
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with the dispersion Uq^- The latter incorporates the effect of a fluctuating magnetic field. 
The 'impurity' part and the effective baths are mixed by i/mix; where the electron (spin) 
of the 'impurity' couples to a fermionic (bosonic) bath via V (/). The parameters for 
the effective baths, i^^o-, 0;^^, V and / are determined self-consistently by equating the 
Green's function and the spin susceptibility obtained in the effective impurity model 
with those in the original lattice model. The local Green's function for c electrons 
and the susceptibility for / moments are defined as Geo- (t) = — (^TT-Cjo-(r)c|^(0)^ and 

Xfa (t) = (TtS{^{t)S-^{0)^ , which are calculated by solving the impurity model. We 
solve the effective Hamiltonian i^EDMFT by means of NCA [32-34] and calculate G*^^ and 
xj^^ [35]. The Green's function and the susceptibility satisfy the following equations, 

1 -1 



(6) 



X}7M = [GU^) + M^iu)]-\ (7) 
where Gaa{^) and Gha{uj) describe the effective fermionic and bosonic bath: 

(9) 



Gha{'^) = ^ 



iU — iUqa U) + UJqa 



9 

We can obtain the momentum-independent self-energy and Ma with these equations. 

The local Green's function and susceptibility in the lattice system are obtained in 
terms of the self-energy Eg. and by using the relations 

Gt\^) = I de -—- (10) 

where Pf(e) and pj{e) are the density of states for the free c-electrons and bosons 
representing spin fluctuations, which are defined as pt{£) = jj^k^i^ ~ ^k) and 
pj{e) = jj: "^qS^e — Jq). Equations (6)-(ll) complete the self-consistent loop. 

3. Results 

We first discuss how the RKKY-exchange among /-spins renormalizes the one-particle 
spectrum and the spin dynamics at i? = 0. One-particle Green's function Gca{k,uj) for 
conduction electrons is calculated with the self-energy So-(c<j), 

GUk, 00) = ^ ^1 (12) 

uj + fi+ ^aB - Sk - 2^a{uJ) 

In Fig. m we plot the one-particle spectral function Ac{k,uj) = —lmGca{k,uj)/ir for 
different values of the RKKY-interaction J at a given temperature T = 0.05. In the 
absence of J, the spectrum acquires an insulating gap resulting from the Kondo coupling 
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Figure 1. The momentum 
resolved one-particle spectra 
Ac{k,uj) at Jk = 0.5, T = 0.05 
for different values of the RKKY- 
interaction J. The data are 
plotted along the direction 
in the fc-space. 



Figure 2. The dynamical spin 
susceptibility, Imx/((7,tj), at Jk = 
0.5 and T = 0.05 for different 
values of the RKKY-interaction J . 



Jk, since we deal with the half-filled band here. The size of the gap A can be defined 
as the energy difference between two peaks with the lowest excitation energy, which are 
located at = [0, 0, 0] and k = [n, n, n]. The gap is almost unchanged upon introducing 
the RKKY-interaction. This result partly comes from the NCA approximation, which 
has a tendency to overestimate the strength of Kondo singlets. We see, however, that 
the sharp peaks located around k = [7r/2, 7r/2, 7r/2] are suppressed and shifted to lower 
frequencies as the RKKY-interaction increases, signaling that the local singlet state 
gradually gets unstable in the presence of the RKKY-interaction. 

The effects of the RKKY-interaction are observed more clearly in spin excitations. 
We calculate the dynamical spin susceptibility XfaiQ,^^) with two-particle self-energy 

Xfaiq,uj) = —j^——. (13) 

In Fig. 121 we plot the imaginary part of the dynamical susceptibility lmxf{q,uj) for 
different values of J. The dominant structure in Imx/(g, cj) is the low- frequency peak 
featuring spin-triplet excitations. As the RKKY-interaction increases, the spin-triplet 
peak at g = [vr, vr, tt] is strongly enhanced and the position of the peak shifts to 
lower frequencies, indicating that the system gradually approaches the antiferromagnetic 
instability. 
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We now discuss how the effects of the RKKY-interaction are modified in the 
presence of a magnetic field. In Fig. El the one-particle spectra for up-spin electrons, 
Ac]{k, ui), are shown aX B = 0.3. We recall that the Zeeman splitting due to the magnetic 
field pushes Ac^ik^u) down to the negative- frequency side, so that the peak position at 
k = [0, 0, 0] is shifted below the Fermi level, and thus the insulating gap is closed. This is 
indeed observed in the results for J = in Fig. El Upon introduction of J, there appear 
two notable features. First, the effective Zeeman shift is increased in the presence of J, 
which is seen in the high frequency region of Fig. (HI This is caused by the competition 
between the Kondo and RKKY couplings. For J = 0, the Kondo-singlet formation 
suppresses the effect of the Zeeman shift. Since the RKKY interaction enhances spin 
fluctuations of /-electrons, it makes the Kondo effect weaker, thus reviving the large 
Zeeman shift. 

Another important feature is the characteristic low-energy behavior. For example, 
in the low-energy region for J = 0.3, the dispersion is highly renormalized to form a fiat 
band around ~ 0, which cannot be explained in terms of the effective Zeeman shift 
mentioned above. Such behavior is seen more clearly in the quasi-particle dispersion of 
c-electrons, LJca{k), plotted in Fig. jH This remarkable effect is regarded as the mass 
enhancement caused by critical /-spin fluctuations. In fact, the system at J = 0.3 is 



B = 0.0 

B=Q.3 up-spin 

B=0.3 down-spin 




Figure 3. Tlie one-particle spec- 
tral function for up-spin (majority- 
spin) electrons, Ac^{k,uj), at the 
magnetic field B — 0.3: Jk — 0-5, 
T = 0.05. 



Figure 4. The quasi-particle 
dispersion ujccr{k) for (a) J = 0.0 
and (b) J = 0.3. The dispersion 
is obtained by tracing the peak- 
position of the one-particle spectra 
Acik,uj). 
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/=0.0 




Figure 5. The magnetic field dependence of the gap for different values of J. 

very close to the quantum phase transition, as discussed momentarily below. Conduction 
electrons scattered by such critical /-spin fluctuations get heavier, resulting in the fiat 
dispersion. The corresponding gap is shown in Fig. El as a function of the field. In the 
absence of the RKKY-interaction, the gap is reduced linearly and closed at the field 
corresponding to the gap size B ~ 0.16. As the RKKY-interaction is increased, the gap 
can persist up to larger fields reflecting the renormalization effect mentioned above. 

In order to see how /-spin fluctuations behave in magnetic fields, we plot the field 
dependence of the imaginary part of the dynamical susceptibility lmxfa{,<i = [tt, vr, vr], cj) 
(a = x,z) in Figs. IHl and Comparing the susceptibility for J = 0.1 and 0.3 at 
B = 0, one can see that [tt, tt, tt] triplet modes are softened in the presence of the 
RKKY interaction. As the magnetic field is introduced, the longitudinal (2;-component) 
and the transverse (x-component) susceptibilities exhibit different properties. Since 
Imx/z(g = [tt, TT, tt], cij) excites a.n Sz = triplet mode, its peak position is mostly 
unchanged even at finite fields, as is seen both for the J = 0.1 and 0.3 cases. On the 
other hand, the transverse susceptibility Imx fx{(l = [tt, tt, tt],uj) shows a splitting in the 
spectrum, corresponding to Sz = ±1 triplet modes. The peak in the higher-energy 
side is smeared considerably by the life-time effect. It should be noticed here that 
the peak structure, which is shifted to lower frequencies, gets sharper as B increases, 
and this tendency becomes more significant as J becomes large. For example, for the 
J = 0.3 case, although weak magnetic fields {B = 0.1) merely suppress spin fluctuations 
and thus reduce the spin-triplet peak, larger magnetic fields {B = 0.2) increase 
the peak intensity again, signaling that the system approaches the antiferromagnetic 
instability. To see the antiferromagnetic instability directly, we show the transverse 
static susceptibility of localized /-spins, Xfxil = [tt, tt,tt]), in Fig. |H1 For smaller RKKY 
couplings (J = 0.1,0.2) there is no indication of the instability as the magnetic field 
increases. We can see the rapid increase of XfxiQ = [7r,7r,7r]) for J = 0.3 reflecting the 
antiferromagnetic instability: the system is driven to the transverse antiferromagnetic 
phase {x-y plane) with finite values of z-components (canted order). 
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Figure 6. The imaginary part 
of the dynamical susceptibility 
Imx/a(9 = [7r,7r,7r],tj): (a) z 
direction and (b) x direction for 
different magnetic fields at J = 0.1. 



Figure 7. The same plots of the 
dynamical susceptibility as in Fig. 
El The RKKY interaction is chosen 
as J = 0.3. 




Figure 8. The transverse spin susceptibility in the x direction Xfxil = [tt, tt, tt]) as a 
function of the magnetic field B for different strength of the RKKY-intcraction J. 



The above characteristic properties also affect the static quantities in magnetic 
fields. For instance, we plot the magnetization curve in Fig. IHl At J = 0, Mf (Mc) has 
a relatively large (small) response to the field. As the strength of the RKKY interaction 
increases, Mf is suppressed whereas Mc is enhanced. Although the sharp spin-gap 
behavior in the magnetization is somehow obscured for the parameters exploited here, 
the above magnetic-field dependence indeed comes from the competition between the 
Kondo coupling and the RKKY interaction, as explained below. The increase of J has 
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J^O.O / = 0.1 



/ = 0.2 



/ = 0.3 




Figure 9. The magnetization process: (a) c-electrons, (b) /-electrons and (c) the 
total contribution. Mc, Mf and Mtot are defined by Mc — {S^), Mc — {Sf), and 
Mtot = Mc + Mf. 




Figure 10. The density of states for up-spin conduction electrons for some different 
choices of the magnetic field B and the RKKY-interaction J. 



a tendency to suppress the Kondo singlet formation, so that the conduction electrons 
start to polarize rather freely for finite J, increasing the magnetization. On the other 
hand, Mf is generally suppressed with the increase of J. Such behavior is consistent 
with the change of the density of states (DOS) shown in Fig. ^1 Comparing the DOS 
between the cases of J = and 0.3, we see that the effective Zeeman shift in conduction 
electrons is somewhat enhanced for J = 0.3, which is clearly observed around ~ 1. 
This shift gives rise to the increment of Mc as seen in Fig. IHl 

Finally, we make a comment on the low-energy behavior of the DOS in Fig. ^| If 
we look at the DOS around ~ 0, the shift of the peak structure is smaller for J = 0.3 
than J = contrary to the effective Zeeman shift. This is what we have claimed as the 
mass renormalization effect due to critical transverse /-spin fluctuations in Fig. |H1 it is 
a typical correlation effect between conduction electrons and /-spins, which is beyond 
the naive mean-field picture. 
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4. Summary 

We have studied the KLM with the RKKY exchange interaction among local /-spins 
in the presence of a magnetic field. For this purpose we have exploited an extended 
version of DMFT combined with NCA. This approach enables us to treat inter-site spin 
correlations dropped in the DMFT framework. By investigating the one-particle spectral 
function and the dynamical spin susceptibility we have confirmed that the RKKY- 
interaction has a tendency to enhance field-induced magnetic instability to the transverse 
magnetic order of /-spins, as should be expected. We have clarified two characteristic 
features caused by the competition of the Kondo and RKKY interactions in magnetic 
fields. One is a simple mean-field type effect: the RKKY interaction enhances the 
Zeeman shift for conduction electrons, which modifies the profile of the one-particle 
spectrum at high frequencies as well as the magnetization curve. Another important 
feature is a renormalization effect observed in the low-energy region: the quasi-particle 
spectrum gets renormalized to form heavy quasi-particles when the system is close to 
the phase transition point. This effect appears as if the RKKY interaction protects 
the Kondo insulating gap in magnetic fields. Anyway, this remarkable feature at low 
energies is a typical correlation effect beyond the mean field picture, which characterizes 
the critical region of the field-induced phase transition in the Kondo insulator. 

Field-induced transitions to the transverse antiferromagnetic order should be 
observed in heavy-fermion Kondo insulators. Unfortunately, there have been few 
examples of such transitions reported so far. Recently, it has been experimentally 
found [49, 50] that a typical fiUed-skutterdite, CeOs4Sbi2, which may be classified as 
an Kondo insulator, exhibits a field-induced magnetic phase transition. This transition 
could be a transition to the transverse antiferromagnetic phase discussed here, although 
more detailed study should be done since the above skutterdite has rather complicated 
electronic structure including orbital degrees of freedom. We think that a variety of 
compounds exhibiting such field-induced transitions may be found in the class of heavy 
fermion Kondo insulators. 
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